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Possible explanations of the observed accelerated expansion of the Universe are the introduction 
of a dark energy component or the modifications of gravity at large distances. A particular difference 
between these scenarios is the dynamics of the growth of structures. The redshift distribution of 
galaxy clusters will probe this growth of structures with large precision. Here we will investigate 
how proposed galaxy cluster surveys will allow one to distinguish the modified gravity scenarios from 
dark energy models. We find that cluster counts can distinguish the Dvali-Gabadadze-Porrati model 
from a dark energy model, which has the same background evolution, as long as the amplitude of 
the primordial power spectrum is constrained by a CMB experiment like Planck. In order to achieve 
this, only a couple of hundred clusters in bins of width = 0.1 are required. This should be easily 
achievable with forthcoming Sunyaev-Zel'dovich cluster counts, such as the South Pole Telescope in 
conjunction with the Dark Energy Survey. 



I. INTRODUCTION 

During the course of the last ten years evidence for an 
accelerated expansion of the Universe has been mount- 
ing. Initially driven by Type la Supernovae (SNe) obser- 
vations by the Supernovae Cosmology Project (SCP) and 
the High-z Supernovae Search Team [1-3], the combina- 
tion of large scale structure surveys like the 2 degree field 
(2dF) and the Sloan Digital Sky Survey (SDSS) [4, 5], 
cosmic microwave background (CMB) [6, 7] and x-ray 
observations of clusters of galaxies [8, 9] confirmed these 
findings. In recent years even larger samples of Type 
la SNe [10-12] indicate strongly that the expansion of 
the Universe is speeding up. While the observations are 
getting better and better, theoretical models of acceler- 
ated expansions are still in its infancy. The most natural 
candidate seems to be a cosmological constant, hovifever 
this comes with a drawback of fine tuning of the initial 
conditions to unnatural values. 

In general there are three possible approaches to 
obtain accelerated expansion of the Universe. Firstly 
there could be an additional component in the energy- 
momentum tensor, which dominates the late time 
dynamics of the Universe and leads to accelerated 
expansion. This component requires negative pressure 
and an example of such a fluid could be a nearly homo- 
geneous scalar field prevailing in the Universe [13-17]. 
The second possibility, is that Einstein's equations of 
gravity require modification on large scales [18-22]. 
The third possible explanation is that the Universe is 
inhomogeneous on large scale, although it is hard to 
explain the boundary conditions for such a scenario in a 
natural way [23]. In this article we concentrate on the 
first two possibilities. 

Numerous dark energy models have been suggested in 
the past, where canonical scalar fields with various poten- 
tials are only a sub-class [13-17]. In all these models the 
background evolution can be described by the equation 
of state, i.e. the relation between the pressure and den- 
sity of the dark energy fluid. The standard dark energy 



models also have a standard evolution of the growth of 
structures. Theories of modified gravity are only emerg- 
ing as possible explanations for the observed acceleration. 
One of the most prominent models is one, which is moti- 
vated by brane cosmology involving extra dimensions and 
hence lead to an effectively modified Friedman equation 
on the 3-1-1 dimensional brane [18, 21]. Another possibil- 
ity is to have inverse curvature terms in the Lagrangian 
of gravity. The additional terms become relevant at the 
late time evolution of the Universe and can lead to ac- 
celerated expansion and hence explain the Supernovae 
data [19, 20, 24]. When modifying gravity, one has to be 
very careful so as not to violate high precision tests of 
gravity in the local Universe, like the Solar System, and 
also not to introduce unphysical features in the theory, 
like negative energy eigenstates, namely ghosts, or super- 
luminal modes [21, 25-32]. For the rest of this article we 
will concentrate on the Dvali-Gabadadze-Porrati (DGP) 
model as an example for a modified gravity model. The 
main reason for this is that in this model the dynamics 
of linear perturbations has been worked out [33-35] . Al- 
though, some care is required on large scales [34, 36] , this 
is not important for the cluster count test, which probes 
power spectra on sub-horizon scales. 



II. DARK ENERGY VERSUS DGP 
COSMOLOGY 

For the analysis presented here we take a DGP model 
as our fiducial model. The DGP model is an example of 
a brane world model, where our four dimensional Uni- 
verse resides in a five dimensional space on a brane [18]. 
The novelty about this model, compared to other brane 
world models, is that the flat extra dimension is allowed 
to be large [21]. The changes of the cosmological evolu- 
tion on the four dimensional brane, compared to standard 
gravity, are caused by 'leakage' of gravitational degrees 
of freedom off the brane. Cosmological solutions allow 
empty Universes with accelerated expansion [37, 38] . The 
resulting modifled Friedman equation for the simplest of 
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these models is given by [21] 
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where i/oGP is the Hubble parameter in the DGP model, 
G the four dimensional gravitational constant and p the 
energy density of the constituents of the Universe. For 
the analysis in this article we focus on the late time evolu- 
tion and set overall density p equal to the matter density 
Pm- Gravity on the brane is affected by 5-dimensional ef- 
fects above the crossover scale r-j,. Self-acceleration of the 
Universe takes place today if rc ^ Hq^. For a flat Uni- 
verse, the parameter rc is related to the matter contents 
n^, by [21] 
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which ensures that one obtains Hq as the expansion rate 
of the Universe today. 

In order to compare the evolution of Universe in the 
DGP world with the one of a dark energy model, we 
note that the Hubble parameter in standard gravity at 
late times is given by 
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with w(a) the equation of state factor of the dark en- 
ergy component. By identifying the Hubble parameters 
in Eqn. (1) and (3) we obtain w{a) which mimics the 
background evolution of the DGP model 



w{a) = -1 + 
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(4) 



with J] = y/n„ia^'^ + l/(2rciJo)^. This tends to w{a = 
1) = -1/(1 + f^m) today and to w{a 0) = -1/2 at 
early times, if we include the condition in Eqn. (2). In 
Fig. 1 we show the mimicking equation of state factor 
(dot-dashed line), if we set fim — 0.3. Given that a 
ACDM model {w — —1, dashed line) fits current observa- 
tions very well, this model might be already under some 
strain from observations. The behaviour of the mimic 
model at late and early times also allows to obtain the 
coefficients of a popular parameterization of the equation 
of state factor, which is given by w{a) = wq + Wa{l — a) 
[39, 40]. For the case we discuss here the coefRcients 
are: Wq — —0.77 and Wa = 0.27. The equation of state 
factor for this parameterization is shown by the triple 
dot-dashed line in Fig. 1. 

We will now address the question how a purely ge- 
ometrical probe, like the Supernovae magnitude - red- 
shift relation, can constrain the DGP model. For this 
we calculate the model predictions for the magnitudes, 
which are given by m oc logdi, and the luminosity dis- 
tance cIl — {I + z) dz' /H{z'). In Fig. 2 we plot the 
magnitude difference to the fiducial DGP model with 
Hq = 72km/s/Mpc and fim = 0.3. From this we obtain 
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FIG. 1: The mimicking w{z) for the DGP model (dot-dashed 
line). The triple dot-dashed line is the 'best fit' parameteri- 
zation with 'w{z) — wo + Waz/{1 + z). We find wo = —0.77 
and Wa ~ 0.27 for the fiducial DGP model. 
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FIG. 2: Magnitude-redshift relation with respect to the fidu- 
cial DGP model with rc — lAS/Ho (solid). The shaded region 
is for models where rc varies from 1.38/Ho to lAS/Ho- The 
errorbars are the statistical errors from a SNAP-like experi- 
ment. The dashed line is for ACDM Universe. 



with Eqn. 2 a fiducial value for TcHq 1.43. The shaded 
region shows the magnitude difference if we vary TcHq 
between 1.38 and 1.48. Note that for these values the 
Hubble parameter today is still well within the errorbars 
of its measured value of Hq = (72 ± 8)km/s/Mpc [41]. 
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We also include mock data points from a SupcrNovae 
Acceleration Probe (SNAP)-like survey [42]. Wo assume 
a distribution of 2000 Type la SNe between z = 0.3 
and z = 1.7 with a magnitude error of am = 0.15 for 
an individual SNe [43]. We plot the statistical error- 
bars in redshift bins of Az = 0.1. A SNAP like exper- 
iment will not be able to distinguish DGP models with 
1.38 < rcHo < 1.48. Realistically the situation is worse, 
since such large bins probably hit the systematic error 
limit of the survey. Nevertheless, it is evident from the 
plot, that SNAP could distinguish the DGP models from 
a ACDM model (dashed line) . Of course the dark energy 
model with an equation of state given by Eqn. (4) , would 
be spot on top of the fiducial solid line in the plot. A pure 
geometric probe, like the SNe magnitude-redshift relation, 
can not distinguish a modified gravity model from an ar- 
bitrary dark energy model. However, a SNAP-like survey 
is likely to probe the geometry and hence the background 
evolution very well. 

So far we have concentrated on SNe as a cosmological 
probe and established that SNe alone cannot distinguish 
dark energy from the DGP model. However, recent stud- 
ies have shown that weak lensing, baryon oscillations or 
observations of the integrated Sachs- Wolfe effect have po- 
tential to distinguish dark energy from modified gravity 
[33, 44-47]. The ability of these probes to distinguish 
the DGP model from dark energy is driven by the dif- 
ferent dynamics of the growth of structure. The linear 
density perturbations in the DGP model are governed by 
[21, 33, 35, 46] 
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with /3 ~ 1 — ^HrcO [3/a + H' /H], and primes denoting 
derivatives with respect to the scale factor a. Note, that 
if we set the term 1/(3/3) = 0, we obtain the standard 
gravity result. For the evolution of linear perturbations 
this is the key difference between the DGP model and 
the mimic dark energy model. 



III. GALAXY CLUSTER COUNTS 

Counts of galaxy clusters are very sensitive to the evo- 
lution of the linear growth factor. Clusters are thought 
to be formed in the peaks of the linear density field, if 
the density is above a threshold density of dc = 1.686, 
for a flat matter dominated Universe. This is the basis of 
the Press-Schechter approach for the distribution of clus- 
ters in mass and redshift [48]. The modern approach is 
to measure this function with N-body simulations, which 
result in [49] 

^iz,M) = -0.316%^^ J- 
dM ^ ' M dM GM 

X exp {-[0.67 - log [5{z)ctm] I^'^^} j(6) 

where n is the number density of clusters of mass M, <tm 
is the density fluctuation of a shell with mass M in terms 



of the mean background mass density and 8{z) is the 
growth factor as a function of redshift. It is clear from 
Eqn. 6 that the mass function is strongly dependent on 
the growth factor 5. Although, the mass function might 
be slightly different in modified gravity scenarios [50], 
the dominant contribution is the different dynamics of 
the growth of structures. We will investigate if counts of 
galaxy clusters can distinguish the DGP model from the 
mimic dark energy model. This is similar to the earlier 
works, which concentrate on weak lensing, baryon os- 
cillations and integrated Sachs- Wolfe effect observations 
[33, 44-47]. 

The number of clusters per redshift bin is given by [51] 
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where dV / {dzdO) = [r{z)Y /H is the comoving volume in 
a flat universe, with r{z) = H~^{z')dz' the coordinate 
distance and AO is the angular sky coverage of the sur- 
vey. Afiim(z) is the limiting mass, which will depend in 
general on the parameters of the survey and cosmology. 

The near future will see a few surveys, which will per- 
form a blind search for galaxy clusters and identify them 
via the so called Sunyaev-Zcl'dovich decrement in the 
CMB radiation [52]. The interesting feature of this ef- 
fect, is that it is independent of redshift and hence could 
potentially identify clusters at large redshift. Upcom- 
ing surveys are the South Pole Telescope (SPT) [53], 
Atacama Cosmology Telescope(ACT) [54] and Atacama 
Pathfinder Experiment (APEX) [55]. In order to per- 
form number counts of galaxy clusters in redshift bins 
we also require the redshift. The proposed Dark Energy 
Survey (DES) [56] is designed to obtain the redshifts of 
the clusters identified by SPT. SPT will observe 4.000 
deg^ of sky with a 1 arcminute beam. In order not to 
be too SPT specific in the analysis presented here, we 
do not calculate the limiting mass of the survey from the 
instrument parameters for SPT. We assume an effective 
constant mass limit of = 1.7 x IQ^^h^^MQ, which 
corresponds roughly to the mass limit of the SPT survey 
[51]. However the results we find will apply to any kind 
of galaxy cluster survey. Apart from the cosmological 
parameters mentioned so far, we also require for the ana- 
lytical prediction of the cluster counts the normalization 
and slope of the primordial matter power spectrum. We 
choose a slope of n = 1 and erg = 0.75 [7]. For this survey 
and fiducial model we find « 7200 clusters of galaxies. In 
Fig. 3 we show the results for counts of clusters in red- 
shift bins of Az = 0.1. The solid line is the fiducial DGP 
model, with randomly generated mock data points in- 
cluding Poisson errors. The faint dotted lines on either 
side of the fiducial model correspond to l/{rcHo) = 1.38 
and l/(rc-ffo) = 1-48. Note that these models are hardly 
distinguishable from the fiducial model in the number 
counts. The dot- dashed line is for the mimic dark en- 
ergy model with equation of state given by Eqn. 4. The 
difference is, particularly in the peak of the distribution. 
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FIG. 3: Number of clusters of galaxies above a threshold of 
Miim = 1.7 X 10" /i"^M0, for the concordance ACDM Uni- 
verse (dashed) and the fiducial DGP model (solid). The thin 
dashed line is for a model with w = —0.8. The dotted lines are 
for the Tc = 1.38/-ffo and rc = 1.48/-ffo models respectively. 
The dot-dashed line is for a dark energy model, which mimics 
the background evolution of the DGP scenario. The error- 
bars assume Poisson errors for each bin with Az = 0.1 for a 
SPT-like experiment. The shaded regions are illustrating the 
uncertainties on the mimic model from propagating errors on 
Wo, Wa, and as from SNe and CMB observations. The 
light shaded region is obtained by assuming Sag = 0.03 and 
the dark shaded on for 5 as = 0.01. 



well above the 2a threshold. This is driven by the dif- 
ference in the growth factor of the DGP model, given 
by Eqn. 5 and the mimic dark energy model, which is 
obtained for /3 ^ oo. So prospects are very good that 
galaxy cluster counts can distinguish the modified gravity 
models from dark energy in conjunction with a geomet- 
rical probe like SNAP. However, we also have to address 
the problem of parameter degeneracies. It is very likely 
that the mimic model with a different value of as is de- 
generate with the DGP model. In order to study this 
effect we perform a simplistic error propagation accord- 

ing to SANiz^r = with 
9j — {iuq, Wq, f^ni, cg}. We adopt here the parameteri- 
zation w(a) — Wo + ii'a(l — a) of the mimic model with 
Wo = -0.77 and Wa = 0.27. From a SNAP-Uke survey 
we expect 8wq = 0.05 and bwa = 0.2 with — 0.03 
[43] . Current constraints on a% from combining WMAP3 
and 2dF or SDSS are dag, = 0.03 [7]. This case is illus- 
trated as the light shaded region in Fig. 3. For Planck 
the constraint on the amplitude of the primordial power 
spectrum is likely to improve by a factor 5 compared 
to WMAP [57]. We hence adopt conservatively a prior 
on cTg from Planck, possibly in addition to a large scale 



structure measurement like SDSS, of Jcg = 0.01. This 
is the dark shaded region in Fig. 3. We see from the 
Figure that even if we marginalize over parameter de- 
generacies we can still distinguish the mimic model from 
the DGP model significantly above the la level. One 
problem we have not addressed in this analysis, is the 
uncertainty in the mass limit of the survey [51, 58]. This 
will effectively increase the errorbars in Fig. 3. However, 
self-calibration methods [58, 59] and optical calibration 
of the cluster masses [56, 59, 60] are very promising and 
could constrain the error in the mass limit. In order to 
obtain some insight in how different the mimic dark en- 
ergy model is from the fiducial DGP model, we also plot 
the number count prediction for the ACDM model (thick 
dashed line) versus a model with w = —0.8 (thin dashed 
line). With self-calibration it should be possible to dis- 
tinguish these two models [56, 59], particularly if priors 
from the Planck CMB mission on the other cosmological 
parameters are available [56, 57]. Given that the differ- 
ence between these two models is smaller than the one 
between the DGP and mimic dark energy model, clus- 
ter counts should be able to distinguish modified gravity 
from dark energy. 



IV. CONCLUSION 

We have shown that the background evolution of the 
DGP model can be entirely mimicked by a dark energy 
model with a suitably chosen equation of state, given in 
Eqn. 4. Due to this fact Type la SNe alone can not 
distinguish modified gravity models from dark energy. 
However, due to the different time evolution of the lin- 
ear perturbations, galaxy cluster counts can distinguish 
the DGP model from a dark energy model with the same 
background evolution. The difference between these two 
models is more then at the 2a level, even for relatively 
low number of clusters « 7000. This result is only slightly 
weakened if we include parameter degeneracies in the 
mimic model. We have presented a very simplified anal- 
ysis in the article, but the striking difference between 
the mimic dark energy and the DGP model, raises hopes 
that this result also stands up in a more careful analy- 
sis, which for example includes uncertainties in the gas 
physics of clusters in order to establish the mass limit 
[61]. 

A further interesting question would be to see whether 
our findings also hold up for other modified gravity sce- 
narios, like the inverse curvature model [19]. In order 
to answer this question a perturbative analysis of these 
models is required, which is not an easy task. It would be 
very surprising if the standard perturbation result holds 
for these models. In general we would expect that mod- 
ified gravity models lead to a modification of the growth 
of structure. Clusters of galaxies in combination with 
a geometrical probe, like SNe distances, could provide 
a uniquely sensitive window for the distinction between 
general modifications of gravity and dark energy. 
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